Self- Similar Hot Accretion Flow onto a 
Rotating Neutron Star: Structure and 

Stability 

Mikhail V. Medvedev 1 and Ramesh Narayan 2 

1 CITA, University of Toronto, Toronto, Ontario, M5S 3H8, Canada 
2 Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138 



Abstract. We present analytical and numerical solutions which describe a hot, vis- 
cous, two-temperature accretion flow onto a rotating neutron star or any other rotating 
compact star with a surface. We assume Coulomb coupling between the protons and 
electrons, and free-free cooling from the electrons. Outside a thin boundary layer, 
where the accretion flow meets the star, we show that there is an extended settling 
region which is well-described by two self-similar solutions: (i) a two-temperature so- 
lution which is valid in an inner zone r < 10 2 ' 5 (r is in Schwarzchild units), and (ii) 
a one-temperature solution at larger radii. In both zones, p oc r~ 2 , Q oc r~ 3 / 2 , v oc 
r° 5 T p oc r" 1 ; in the two-temperature zone, T e oc r -1 / 2 . The luminosity of the settling 
zone arises from the rotational energy of the star as the star is braked by viscosity. 
Hence the luminosity and the flow parameters (density, temperature, angular velocity) 
are independent of M. The settling solution described here is not advection-dominated, 
and is thus different from the self-similar ADAF found around black holes. When the 
spin of the star is small enough, however, the present solution transforms smoothly to 
a (settling) ADAF. 

We carried out a stability analysis of the settling flow. The flow is convectively 
and viscously stable and is unlikely to have strong winds or outflows. Unlike another 
cooling-dominated system — the SLE disk, — the settling flow is thermally stable 
provided that thermal conduction is taken into account. This strong saturated-like 
thermoconduction does not change the structure of the flow. 



THE SETTLING FLOW 

At small mass accretion rates, < 10~ 2 of the Eddington rate, black holes (BHs) 
accrete via an ADAF — a hot, two-temperature, radiatively inefficient, geometri- 
cally thick, advection-dominated accretion flow [1,2]. In contrast, accretion onto 
compact stars, e.g., a neutron star (NS) may occur via either an ADAF of a settling 
solution [3]. The latter corresponds to strongly rotating stars only. In the settling 
flow the rotational energy is extracted from the star via viscous torques in the 
boundary layer where the flow meets the star surface. The extracted energy heats 



the flow and ultimately escapes from the flow as free-free radiation. In addition, 
viscosity extracts angular momentum from the star as well. In the stationary flow, 
this angular momentum must be transported though the flow to the outermost 
radii, where it goes into the ambient medium. It is this huge angular momen- 
tum flux J which modifies the entire structure of the accretion flow and makes 
it drastically different from an ADAF. Note that no angular momentum may be 
extracted from a BH horizon by viscosity. Thus a viscous settling flow may exist 
in NS systems and not in BH systems. 

The structure of the steady, rotating, axisymmetric, quasi-spherical, two- 
temperature settling flow has been found analytically and confirmed numerically [3] . 
We use the height-integrated form of the viscous hydrodynamic equations with the 
Shakura-Sunyaev-type viscosity parametrized by dimensionless a. We assume vis- 
cous heating of protons, Bremsstrahlung cooling of electrons and Coulomb energy 
transfer from the protons to the electrons; we neglect Comptonization but include 
thermal conduction in the form discussed in the next section. In the inner zone 
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FIGURE 1. Profiles of density p (g cm~ 3 ), proton temperature T p (°K), electron temperature 
T e (°K), angular velocity SI (in units of the Keplerian angular velocity at the NS radius Rns), 
and radial velocity v (in units of c) for accretion flows with a = 0.1 and (m, s) = (0.01,0.3) - solid 
line, (0.0001, 0.3) - short-dashed line, (0.01,0.1) - medium-dashed line, (0.01,0.7) - long-dashed 
line. 



r < 10 2 - 5 (r is in Schwarzchild units, R s = 2GM/c 2 ), the flow is two-temperature 
with the density, proton and electron temperatures, angular and radial velocities 
scalings as 

p = p r~ 2 , T p = T pQ r'\ T e — T e0 r _1,/2 , Q = Q r~ 3 / 2 , v = v r°, (1) 

where po, T pQ , T e0 ,Q ,vo are functions of M, a and the star spin s = Q*/Qx(R*), 
and Qk(R) = (GM/i? 3 ) 1//2 is the Keplerian angular velocity. In the outer zone r > 
10 2 ' 5 , we have T e = T p oc r~ l and the same other scalings. This self-similar solution 
is valid for the part of the flow below the radius r s related to the mass accretion 
rate m (in Eddington units, M Edd = 1.4 x 10 18 m g/s, and here m = M/Mq): 

m < 2.2 x l(T 3 a 2 lS 2 3 r- 3 1/2 , (2) 

where r s 3 = r s /10 3 , a .i = a/0.1, etc.. The numerical solution of the hydrodynamic 
equations with appropriate inner and outer boundary conditions is represented in 
Figure 1. It is in excellent agreement with the self-similar solition (1). 

There is a remarkable property of the settling flow: all quantities, e.g., p, T p , etc., 
except v, are independent of the mass accretion rate M. This happens because the 
gravitational energy of the accreting gas is much smaller than the energy extracted 
from the rotating star. It is this energy which dominates the flow luminosity. We 
should also remark that the settling flow is more similar to a steady, radiative 
cooling-dominated "atmosphere" rather than to a rapidly infalling flow: the radial 
infall velocity is constant and is much less then the free-fall velocity v/vg oc r 1 / 2 — > 
as r — > 0. The structure of the settling is very sensitive to the rotation rate of the 
central star. As the angular velocity decreases below few percents of the Keplerian 
value, the settling flow smoothly transforms into a conventional ADAF solution, as 
represented in Figure 2. 

It was shown that the settling flow is (i) convectively stable and (ii) may not 
have strong winds and outflows (the Bernoulli number is negative) if the adiabatic 
index satisfies 

3(1 -s 2 /2) 
7 > — — ~ 1.5 

7 (2-s 2 /2) 

Other properties of the settling flow are discussed elsewhere [3,4]. 

THERMAL STABILITY OF THE SETTLING FLOW 

The settling flow is cooling-dominated. Thus it is similar to the Shapiro- 
Lightman-Eardley solution [5] which is known to be thermally unstable [6] . Hence, 
our settling solution may be unstable as well. To study the thermal instability in 
sheared circular accretion flows we use the shearing sheet approximation with the 
velocity given by 



(3) 



V (a;) = 2Axy, 



(4) 



where 2 A = dV /dx is the shear frequency, x, y, z are the radial, azimuthal, 
and vertical coordinates, and "hat" denotes a unit vector. We assume that there 
is a Coriolis acceleration due to ft — z. The vorticity and epicyclic frequency 
are then 2B = 2A + 2fi, K% pi = AQ.B. For a Keplerian-type flow, Q, oc iT 3/2 , 
which is the case for both the settling and SLE solutions considered below, one 
has 2A = —(3/2)1) and 2B = 12/2. We assume that perturbations have only s- 
component, which corresponds to axisymmetric perturbations. We ignore motions 
in z direction. We use hydrodynamic equations with thermoconductive flux. 

The settling flow is hot (sub-virial), so that the mean-free-path (of both electrons 
and protons) is larger than the size of the system. Hence the conventional Spitzer 
theory fails. Without collisions but in the presence of magnetic fields electrons 
stream freely along the field lines, therefore the parallel heat flux remains large. 
In contrast, transverse heat flux is greatly reduced in a magnetic field because 
electrons are tied to the field lines on the scale of the Larmor orbit and cannot 
move across the field lines too far. In a tangled field, however, electrons can jump 
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FIGURE 2. Same as in Fig. 1 for 7 = 4/3 and s = 0.3 (solid curve) and s = 0.01 (dotted 
curve). The self-similar slopes for an ADAF flow and a settling flow are shown for comparison. 
The long-dashed curves represent the same solution as the solid curve, but with ten times higher 
temperature at R ou t- 



from one field line to another and thus transfer heat across the field [7]. The average 
thermal conductivity in tangled fields is 

k b ^ nkEVrlB^ — 10~ 2 nkBVTR^-i'&-i, (5) 

where Ib = £-R is the correlation scale of the magnetic fields set by flow turbulence 
and •& < 1 takes in to account that only a fraction of all particles may pass though 
the magnetic barriers. The rest of them will remain trapped in magnetic wells and, 
hence, will not transport energy to large distances Ib- The thermoconductive 
flux is 
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where a r 



W-^-i^-i F(e,p) (6) 



Qk dx 

is the "alpha-prescription" thermo-conductivity coefficient and F(e,p) [1 < 
F(e,p) ^ 15] takes into account that conduction may be dominated by the protons 
or the electrons, depending on the flow conditions. Here we used that vt — c se and 
H/R ~ c s /vfi ~ c s /VL K R, where H ~ R is the accretion disk scale height. 

Now, from hydrodynamic equations for the perturbations, we obtain the disper- 
sion relation 
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In the large-/c limit, (7) yields the growth rate of the thermal mode. The stability 
criterion Imu < is (for as 2 < a c ) 
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that is, thermal modes with kR > 2 are stable. Whether the mode kR = 1 is stable 
or not cannot be reliably determined from our local approach. A global stability 
analysis is necessary to properly account for the effects of geometry and curvature 
on the eigenmode structure. Finally, we conclude that the settling flow is very 
likely stable. 
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